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The alignment of vorticity and scalar gradient with the eigcndirections of the rate of strain tensor is 
investigated in turbulent buoyant nonpremixed horizontal and vertical flames. The uniqueness of a 
buoyant nonpremixed flame is that it contains regions with distinct alignment characteristics. The 
strain-enstrophy angle 'F is used to identify these regions. Examination of the vorticity field and the 
vorticity production in these different regions indicates that 'T and consequently the alignment 
properties near the flame surface identified by the mixture fraction band F~F>, differ from those in 
the fuel region, F>F SI and the oxidizer region, F<F sl . The F=F J( band shows strain-dominance 
resulting in vorticity/or alignment while F>F SI (and F<F SI for the vertical flame) band(s) show(s) 
vorticity//? alignment. The implication of this result is that the scalar dissipation, e F , attains its 
maximum value always near F^F st . These results are also discussed within the framework of 
recent dynamical results [Galanti et al., Nonlinearity 10, 1675 (1997)] suggesting that the 
Navier-Stokes equations evolve towards an attracting solution. It is shown that the properties of 
such an attracting solution are also consistent with our results of buoyant turbulent nonpremixed 
flames. © 1998 American Institute of Physics. [S1070-6631(98)01709-7] 
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I. INTRODUCTION 

More than 40 years ago, Batchelor and Townsend 1 sug- 
gested that the material lines consisting of fluid particles tend 
to align along the direction of the largest principal rate of 
strain, denoted by a. They also indicated 1 that this alignment 
is unlikely to be perfect since the principal axes of strain 
rotate relative to the fluid. Indeed, Kerr, 2 and Ashurst et al ,, 3 
observed that in homogeneous isotropic and homogeneous 
shear turbulence, the vorticity alignment is not with the larg- 
est strain direction but with the intermediate strain, denoted 
by /?. Dresselhaus and Tabor 4 showed analytically that the 
competition between the strain and rotation determines 
whether the material (or vorticity) lines will align with a or 
/3 direction. Nomura and Elghobashi, 5 and Boratav et al. 6 
showed that for the variable density flows of nonpremixed 
flames with heat release, the vorticity tends to align along the 
a strain direction. 

The analyses of Dresselhaus and Tabor 4 and Boratav 
et al 6 led to a vector equation which can be written for the 
most general case as: 

a — D 0 0 \ 

0 /3-D 0 J A- £A 

0 0 y— D ) 

-(il' + C— ft)xA. (1) 

The components of the unit vector A are the direction co- 
sines of the unit material or vorticity element with respect to 
the strain orthonormal eigendirections e a , eg and e r , where 
a, (3, y denote the eigenvalues of the rate of strain tensor 
with the conventional ordering, a>/3>y. The direction of 
the unit vector A is along the material or vorticity element. 


XI' is the rotational velocity vector of the strain basis axes. C 
accounts for the coriolis effect of the baroclinic torque. For 
incompressible flows, the velocity divergence D, (D = a 
+ (3 + y ) , and the baroclinic term C are zero. For the vortic- 
ity element analysis, the cross product XixA is zero. 

Dresselhaus and Tabor 4 examined the material element 
alignment for an incompressible flow and used the notation A 
for A. Boratav et al. 6 examined the vorticity element align- 
ment and used the notation il for the variable A. £ in Eq. (1) 
is a quadratic nonlinear term which contains the eigenvalues 
a, f3 and y. In an incompressible flow, £ is defined as £ 
= £i ine = a\] + f3\l+ y\y For the variable density vorticity 
alignment, our previous result 6 and further analysis show that 
£^ £ uoni city= + pFl\+ yFl\ — D . It can be shown that 
D appearing in the £ definition cancels with that in the first 
term on the right hand side of Eq. (1). Thus, the divergence 
D does not affect A directly, but D affects the alignment 
indirectly via its presence in the vorticity equation. In turbu- 
lent nonbuoyant flames, the density reduction causes D to be 
positive and thus creates a sink of vorticity 5,11 in the reaction 
zone. Accordingly, a region of weak vorticity and strong 
strain is created in the reaction zone of turbulent nonbuoyant 
flames. 6 The implication of this result in turbulent buoyant 
flames will be discussed in Sec. III. 

The present work aims at addressing the following two 
points: First , as numerous studies of the vorticity-strain 
alignment indicate certain universality, we would like to in- 
vestigate whether our results for buoyant flames show the 
same trends, and to understand why these universality trends 
(or lack of) exist. 

Second, we would like to study the alignment of the 
scalar gradient vector VF with the strain eigenvectors, and 
establish whether this VF alignment is conditional on the 
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FIG. 1. Initial configuration. Cross-section of the solution domain, (a) Hori- 
zontal flame, (b) Vertical flame. U and W (see the arrows) denote the direc- 
tion of the mean flow in x and z directions, respectively. 

value of the mixture fraction F. The alignment of V F with 
the strain eigendirections has direct relevance to the chemical 
reaction process in turbulent nonpremixed flames since it 
controls the temporal development of the scalar dissipation 
rate and subsequently the progress of the chemical reaction. 

The motivation for studying the alignment of strain, vor- 
ticity and scalar gradient in a buoyant nonpremixed flame is 
the unique feature of this flow of consisting of regions of 
quite distinct alignment characteristics as will be discussed 
later. 

The paper is organized as follows: We present a brief 
description of the flow, and the solution method of the gov- 
erning equations in Sec. II. We discuss the vorticity-strain 
alignment characteristics in Sec. in. We introduce a quantity, 
'P, to measure the relative magnitude of vorticity with re- 
spect to strain in Sec. ID A. We examine the regions in the 
flow field having different 'J / characteristics in Sec. Ill B. We 
examine the vorticity-strain alignment pdfs in these in Sec. 
HI C. In Sec. Ill D, we study the effects of the vorticity- 
strain alignment characteristics on the scalar gradient-strain 
alignment by examining the evolution equation of the scalar 
gradient vector [Eq. (4)]. We discuss the universality of the 
alignment results in Sec. IV. 

II. FLOW DESCRIPTION 

The two flow configurations (horizontal and vertical 
flames) chosen for the study are shown in Fig. 1. This figure 
shows a cross-section of the three-dimensional solution do- 
main. The first configuration describes a horizontal flame for 
which the initial interface between the parallel and uniform 
mean-velocity streams of fuel and oxidant is perpendicular to 
the gravity vector. The second configuration describes a ver- 
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tical flame for which the initial interface between the fuel ‘ 
and oxidant streams is parallel to the gravity vector. For both 
cases, the gravity vector is in negative z direction. 

The flow is subsonic and the domain is unbounded and 
thus the thermodynamic pressure is assumed uniform in 
space and constant in time. In small Mach number turbulent 
flows with density variations arising from chemical energy 
release, the kinetic energy is small in comparison to the ther- 
mal energy. In order to compute such a flow, the full set of 
compressible equations may be employed. However, this has 
the computational disadvantage of a rather severe time step 
limitation in order to resolve the high frequency acoustic 
waves. In subsonic (small Mach number) flows, the time 
scale of the acoustic waves is much smaller than those asso- 
ciated with the convection processes. Since the acoustic fluc- 
tuations do not interact effectively with the fluid dynamics, 
they can be neglected. Simplifications to the fully compress- 
ible equations can therefore be made based on the small 
Mach number conditions. The resulting governing equations 
are similar to those of McMurtry el al ., 7 except that the 
buoyancy forces are included in our equations. The three- 
dimensional, time dependent, variable density continuity, 
Navier-Stokes and energy equations are solved together with 
the conservation equations of the mass fractions of the fuel 
and oxidizer. The chemical reaction between the fuel and 
oxidizer follows a single-step, irreversible, binary reaction 
with Arrhenius kinetics. The molecular viscosity, mass dif- 
fusion coefficient, thermal conductivity, and the constant- 
pressure specific heat are assumed to be invariant in time and 
space. The boundary conditions are periodic for the xz and 
yz plane boundaries, while a convective outflow boundary 
condition is imposed along the xy plane boundaries. The 
flow field is initialized with a prescribed energy spectrum 
which is proportional to ke~ k where k is the magnitude of 
the wave number vector. 

The governing equations are discretized using a stag- 
gered grid and a semi-implicit second-order finite differenc- 
ing scheme. The source terms in the energy and species 
equations are discretized using the Crank-Nicolson implicit 
method. The Poisson equation for the pressure is solved us- 
ing a FFT combined with a tri-diagonal matrix solver follow- 
ing the algorithm by Schmidt et al. s 

Three different grids with 96 3 , 128 3 and 192X96X96 
points were used for the simulations with an initial Reynolds 
number based on Taylor microscale, R x = 25. Two grids with 
128 3 and 192X96X96 points were used with an initial R\ 

= 35. In all the simulations, the resolution criteria rjk max 
>1.8 is satisfied, where 77 is the Kolmogorov length scale 
and k max is the maximum resolved wave number in the field. 

The value of R x , and the number of grid points were 
prescribed such that the motion at the smallest scales are well 
resolved. This insures that the velocity and scalar gradients 
(including pressure gradients) are well resolved. The resolu- 
tion accuracy is evaluated by comparing the results , of the 
simulations using grids with successive refinement. For the 
nonbuoyant flame, the difference between the values of the 
pressure Hessian [Eq. (8)] obtained from the simulations 
with the two grids: 96 3 and 128 3 , i.e., increasing the resolu- 
tion by a factor of 2.37, resulted in less than 2% change in 
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the values of the pressure Hessian. For the buoyant flame, 
the gradients arc steeper than those in the nonbuoyanl flame. 
In order to resolve these gradients, we placed more points in 
the direction which has the steepest gradients, namely, the 
gravity direction. Here again, the values of the pressure Hes- 
sian from simulations with the two grids: 96 3 and 192 
X96 2 differed by less than 5% near the F sl surface. 

The turbulence is allowed to develop without chemical 
reaction until the velocity derivative skewness reaches a 
value of approximately equal to 0.5. At that time, the chemi- 
cal reaction is allowed to take place between the two non- 
premixed streams of fuel and oxidizer. 

The ranges of dimensionless numbers tested are: 
Damkohler number=1000, 5000, <*; Froude number=7, 10, 
18 and and R x = 25 and 35. This paper will present results 
of only two buoyant flames, one horizontal and the other 
vertical with Damkohler number=5000, and Froude number 
= 1 0. The initial /? x of the horizontal and vertical flames 
equals 25 and 35, respectively. We will refei also to the 
corresponding nonbuoyant flows of the two cases whenever 
necessary. 

All the simulations continued until a non-dimensional 
time / = 6, which equals about three eddy turnover times. The 
simulations are terminated before the expanding flow ap- 
proaches the boundaries and starts to invalidate the imposed 
boundary conditions. 

All the presented results are obtained at r = 5 to insure 
that the maximum values of scalar dissipation, reaction rate, 
and temperature are already attained. 



FIG. 2. Contours of out-of-plane vorticity on a xz plane at the mid-y loca- 
tion, for the horizontal (lame, ( = 5, Fr = 10, Da = 5000. Only part of the 
solution plane is shown in the z (vertical) direction. Black and white corre- 
spond to negative and positive vorticity. The thick black curves correspond 
to a superposition of surfaces where F^F,,. 


'I' = tan~ 1 


Sij^ij 

R ,j R ,j 


(3) 


The definition of 'P in (3) indicates that laTge values of 'P 
(§>45°) are associated with strain-dominated regions, and 
smaller values of 'P (<^45°) denote enstrophy-dominated re- 
gions. Since we are interested in the effects of chemical re- 
action (density variation) on the turbulence structure, we will 
focus our attention on a mixture fraction zone (0.15*£F 
=£0.85) surrounding the stoichiometric reaction surface (F 
= 0.5). 

In order to identify regions with different characteris- 
tics, we first examine the enstrophy characteristics of the 
different regions (Sec. Ill B) and then compute the 'P char- 
acteristics (Sec. Ill C). 


III. RESULTS 
A. Strain-enstrophy state 

The studies of Batchelor and Townsend, 1 Dresselhaus 
and Tabor, 4 and Boratav et al. 6 show that the relative mag- 
nitudes of the strain and rotation terms determine the 
vorticity/strain alignment characteristics. The rotation terms 
[see Eq. (5.4) of Batchelor and Townsend, 1 Eqs. (12) and 
(19) of Dresselhaus and Tabor, 4 Eq. (2) of Boratav et al. 6 or 
Eq. (1) in the present paper] consist of the vorticity, the 
rotation of the strain axes, and for the variable density case, 
the baroclinic vorticity production. For both incompressible 4 
and compressible 6 flows, when the strain is dominant over 
the rotation, the material lines and the vorticity lines align 
along the direction of the maximum strain a. In this section, 
we focus on the relative magnitudes of the strain and vortic- 
ity. The effects of the rotation of the strain coordinates on the 
alignment will be briefly discussed in Sec. Ill C. 

We examine the relative magnitudes of the strain and 
vorticity at each mesh point in a zone of containing the flame 
surface using a “Strain-Enstrophy State ” plane. The ab- 
scissa and ordinate in that plane are the local enstrophy co 
■ oj/ 2= R jjR ij and the mean square strain rate 5,y5, ; . The 
polar coordinates of a given point on that plane are the dis- 
tance, A, from that point to the origin and the counterclock- 
wise angle, denoted by the Strain-Enstrophy angle, '1 / , and 
measured from the abscissa: 9 


B. RijRij characteristics 

In order to determine the enstrophy (/?,y/?, 7 ) characteris- 
tics of the different regions, we examine the vorticity isosur- 
faces, and investigate the importance of different terms in the 
enstrophy equation in these regions. 

We present in Fig. 2 and Fig. 3 the out-of-plane vorticity 
contours for the horizontal and vertical flames, respectively. 
Positive and negative signs are marked to show the counter- 
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FIG. 3. Same as the previous figure, vertical flame. 
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FIG. 4. F-band-avcraged cnstrophy oi-co/2. Curves correspond to: (a),(b) 
Da = 5000, Fr— 10 vertical (dashed) and horizontal (solid) (lames, (c) Da 
= 1000. Fr = 10 vertical (dashed -dotted). 


rotating vortices. The solid lines are the F-isosurfaces sur- 
rounding the F st surface. It is seen from the figures that the 
horizontal flame F sl surface is saddled by quadruples 
whereas the vertical by dipoles. For the horizontal flame, the 
vorticity above the F s , surface is stronger than that below. 
For the vertical flame, the magnitudes of vorticity on both 
sides of the F sl surface are nearly the same. 

Boratav et al. 6 show that the baroclinic torque is the 
main source of vorticity production in buoyant nonpremixed 
flames. This term changes sign across the F sl surface and 
vanishes at that surface because the density gradient changes 
signs across the reaction zone. They 6 also show that the larg- 
est magnitude of fluid velocity is along the gravity direction 
(w, the z-component of the velocity) and occurs at the F sl 
surface (i.e., at the location of minimum density). Away 
from the F sl surface, the fluid density increases, resulting in 
smaller n\ Thus, dw/dz, which is the major contributor to 
V-u changes sign across the F sl surface. For regions of F 
<F i( , is positive, and for F>F SI , it is negative, thus 
resulting in stronger vorticity in the former than in the latter, 
in the horizontal buoyant flame. 

Similar arguments can be made about the vortex stretch- 
ing term, namely, the dominant contributor to the stretching 
term is dw/dz, which changes sign across F sl . Thus the 
vorticity production due to the stretching is mainly in the 
regions of F>F SI , resulting in stronger vorticity there. 

In summary, the strong baroclinic torque creates vortic- 
ity in both regions of F>F st and F<F st . The velocity di- 
vergence and the vortex stretching terms produce stronger 
vorticity in regions of F>F SI compared to that in F<F SI in 
the horizontal flame. There is no such distinction between 
these regions in the vertical flame. The vorticity attains its 
minimum value at the F s , surface, as presented in Fig. 4, 
which shows the F-averaged enstrophy for different values 
of Froude and Damkohler numbers. 

C. 'I' and alignment 

Based on our analysis in the previous section, we clas- 
sify the flow field into three distinct regions corresponding to 
the following mixture fraction values: (i) F — F sl ; (ii) F 
>F SI ; (iii) F<F sl . In order to have sufficiently large 


Probability ' Fr=1 °, oa=sooo 
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FIG. 5. Horizontal flame, / = 5, Fr= 10, Da = 5000. Top: (a) pdfs of the 
strain-enstrophy angle '}' for three /’-bands: Solid: F~ , dashed: F~ and 
dashed-dotted: F°. Bottom: (b) pdf s of the angle between vonicity and the 
eigenvector corresponding to the largest strain a of the rate of strain tensor. 
Same F-bands as above. 


sample size for the statistics, we choose the following three 
F bands to compute the pdf of 'F: (i) 0.45=£F=s0.55, denoted 
as the F° band; (ii) 0.50< FssO.85, denoted as the F + band; 
and (iii) 0.15=£F<0.50, denoted as the F“ band. 

As discussed earlier, in the horizontal buoyant flame, the 
vorticity is small in F°, and large in F + . Also, the vorticity 
in F + is larger than that in F" . Figure 5(a) shows that the 
flow is strain-dominated not only in F°, which is expected, 
but also in F~ which is a manifestation of the fact that vor- 
ticity is small relative to the strain in these regions. On the 
other hand, as seen from Fig. 5(a), for F + , the vorticity 
strength relative to strain increases. 

Figure 5(b) presents the alignment characteristics for the 
same flow of Fig. 5(a). In this figure, the x axis is the cosine 
of the angle between the vorticity and the largest strain di- 
rection, denoted by cos 6 a . The two strain-dominated bands, 
F° and F~ show a alignment trends, the former being stron- 
ger (i.e., larger probability values). On the other hand, the 
F + band in which the vorticity is more dominant than the 
strain, does not show a alignment trends but (3 (this is not 
shown here due to space limitation). 

Figure 6(a) and (b) shows the results for the vertical 
flame, which are very similar to those for the horizontal 
flame, except that the statistics for F~ and F + are nearly 
identical, and thus only F + results are shown. 
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FIG. 6. Vertical flame, / = 5, Fr=10, Da= 5000. Top: (a) pdfs of the 
strain-enstrophy angle 'F for two F- bands; Dashed: F~ . dashed-dotted: F°. 
The other band F~ statistics are very close to those of F~ and therefore arc 
not plotted. Bottom: (b) pdfs of the angle between vorticity and the eigen- 
vector corresponding to the largest strain a of the rate of strain tensor. Same 
F-bands as above. 


It should be mentioned that the strain rotation term ft' is 
not included as part of the rotation terms in the definition of 
'F. Our computations show that this term is small compared 
to the other rotation terms (i.e., vorticity and baroclinic 
term), and thus it is not included in the 'F definition. 

D. Scalar gradient alignment 

Now we discuss how the vorticity alignment influences 
the scalar gradient alignment with the strain eigenvectors. 
We consider the evolution equation for the scalar gradient: 10 

D 1 

— (VF)=-S-VF+- wxVF + 9>, (4) 

where denotes the viscous term. The contribution of the 
second term on the right hand side of Eq. (4) is to move VF 
toward a direction perpendicular to the vorticity vector. And 
as discussed in the previous section, the flow in the F° band 
is strain-dominated and the vorticity aligns with the a eigen- 
vector. Thus, for the F° band, VF will move toward a plane 
containing and e y . In other words, VF will be strained 
only by /3 and y strains. Since |-yj>|/3|, the deformation of 
VF will be mostly along the y direction. In incompressible 
turbulence, simulations by Kerr 2 and Ashurst et al. 3 show 
that the |y|/|/3| ratio is about 4. Our simulations show that the 


Probability 



FIG. 7. Vertical flame, t-5, Fr — 10, Da = 5000. The pdf of the angle 8 a . y 
(see the text for definition). Note that the most likely value is approximately 
equal to 40°. 


average of this ratio in the F° band is much larger than 4. 
Thus, VF is expected to align strongly along e y in the F° 
band. 

On the other hand, for the F + (and F~ for the vertical 
flame) band(s), the vorticity aligns along the /? eigenvector. 
Thus, VF will be on a plane strained only by a and y strains. 
The dominant direction of deformation will depend on the 
relative magnitudes of a and y. We quantify this relative 
magnitude by computing the a-y angle 6 a . y defined as: 
0 a . y =tan _1 |a)/|y(. Figure 7 displays the pdf of this angle for 
the vertical flame. The pdf shows that the most likely angle 
value is around 40°. This is equivalent to a ratio of |y|/|a| 
= 1.19. 

The pdf of the cosine of the angle between V F and e y 
for the vertical buoyant flame is shown in Fig. 8. As dis- 
cussed above, in the F° band, the |y(/|/3| ratio is large, result- 
ing in good alignment between VF and e y as seen in Fig. 8. 
On the other hand, in F + (and F~ for the vertical flame) 
band(s), as seen in Fig. 7, the pdf of the a-y angle has its 
most likely value at an angle of 40° (note that cos 40° =0.77) 
which is consistent with the cosino of die angle between V F 
and y eigendirection given in Fig. 8 in the F + band. 



FIG. 8. Vertical flame, i~5 ,Fr— 10, Do - 5000. pdf s of the angle between 
V F and the eigenvector corresponding to the largest compressive (negative) 
strain y of the rate of strain tensor. F- bands arc: Solid: F°. Dashed: F* . 
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FIG. 9. F-band-averaged scalar dissipation e F for four different runs with 
Da= 1000 and 5000, and Fr = °° and 10. 


In summary, our results show that both in the absence 
and presence- of buoyancy, the F® band in nonpremixed 
flames is strain-dominated, mainly due to the small vorticity 
and large strain in this region. Accordingly, the vorticity (al- 
beit small) aligns along the a direction in F°. Consequently, 
the scalar gradient VF aligns along y in F°. Figure 9 shows 
that for a wide range of simulation parameters (not all are 
shown here), the peak scalar dissipation = V F - V F always 
occurs within F°. 


Fr=10, Da=5000 

Probability (vertical > 



FIG. 10. The angle between the variable-density vortex stretching cr p 
= tij/p. Vu and vorticity £=<i>lp. Only the vertical flame results are shown 
for two different F-bands for brevity. 


a p and £ in the vertical flame for the F° and F + bands; the 
horizontal flame results in these bands are similar. As ex- 
pected, since the vorticity aligns with a in F° and with /3 in 
F + , the cos 0 pdfs have their most-likely values close to 
unity. 

A quantity whose evolution has been examined in 
detail 12,13 is the scalar A s (in Ref. 13, a) defined as: 


IV. UNIVERSALITY OF ALIGNMENT 


We have discussed in the previous sections that if the 
relative magnitudes of strain and vorticity are known, the 
preferential alignment directions of the vorticity and the sca- 
lar gradient vectors can be predicted. However, we did not 
explain why the dynamics moves toward such a state. 

Recently, Gibbon and Heritage 12 and Galanti era /. 13 
suggested that the alignment trends similar to those obtained 
here and in literature 2,3 can be a manifestation of an attract- 
ing fixed point of the Navier-Stokes equations and in this 
_sense._the.alisnment is universal. In this section, we compute 
certain quantities which appear in their 12,13 analysis and 
check whether our results are consistent with the existence of 
such an attracting solution. 

Gibbon and Heritage 12 and Galanti et al. 13 indicated that 
the fixed point in the Navier-Stokes equations is associated 
with the angle 0 between the vorticity and vorticity- 
stretching vector. This angle is given by : 12,13 


(9 = tan 1 


|wXS-ft»j 
to- S- to 


(5) 


When 6=0, the vorticity and stretching vectors are parallel 
and the vorticity is stretched. Also, when the vorticity aligns 
with an eigenvector of the rate of strain tensor, then the 
vortex stretching vector S -to will also align with the vortic- 
ity, resulting in 0=0 or 9=tt. For Burgers vortex, 0 value is 
strictly equal to zero. It is shown 12,13 that cos 9 approaches 
unity as the solutions move to the attracting fixed point. 

The above analysis can be extended to the variable den- 
sity case 12 if £=co/p is used instead of the vorticity to, and 
the variable density stretching, a p =(to/p)-Vu is used in- 
stead of S • to. Figure 10 shows the pdf of the angle between 


A s = 


to-S-to 
to- to 


( 6 ) 


for incompressible flows. For the Burgers vortex, A 5 is equal 
to the applied external strain which is positive. For the at- 
tracting solution given in Refs. 12 and 13, the solution has a 
stable fixed point for A$>0. For the variable density flows, 
this quantity is defined as: 

= O fp)-S-{to/p) 
s ~ (c o/p)-{to/p ) 

The F-averaged A s distributions for the horizontal and 
vertical flames are shown in Fig. 11. We note that A s values 
are positive in all F bands. Also, the largest A s values on the 
average are in the F° band for both the vertical and horizon- 
tal flames, due to the fact that the denominator, {to/p) ■ {to/p) 
in (7) is small compared to the numerator, in that band. 


As = [(co/p) . S . (co/p)] / [ (co/p) . (co/p)] 



FIG. 1 1 . F-averaged A s values for the horizontal and vertical flames. 
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Another scalar quantity of interest 12,13 is A P , which is 
related to the pressure Hessian P. For the variable density 
(lows, the pressure Hessian P and A P arc given by: 12,13 


P =P 


pJj 


d 

dXj 



( 8 ) 


P ~ (to/p)(to/p) 


( 9 ) 


For Burgers vortex, the quantity A P is strictly negative. For 
the attracting solution given by Refs. 12 and 13, as the solu- 
tion moves towards the fixed point, A P becomes negative. 

Our simulations show that Apis negative for all F bands 
for both the vertical and horizontal (lames. A P attains its 
peak negative (minimum) value in the F° band for both 
flames. 

Finally, we examine how buoyancy changes the charac- 
ter of the pressure Hessian tensor P: Eq. (8) shows that for 
constant density, the pressure Hessian tensor P pi j , is a sym- 
metric tensor and its eigenvalues are real. This is not the case 
for the variable density flow. As an example, we write the 
following two elements of P p ij : 


d ■ 

I 1\ 

dP 

l_ 

i 

d 2 p 

dx 

\PI 

dz 

p 

dxdz ' 

d | 

' l ~\ 

dP 

i 

d 2 P 

dz ' 

tP/ 


p 

dzdx ' 


( 10 ) 


and the two elements are not necessarily equal in the pres- 
ence of buoyancy. Therefore, buoyancy creates eigenvalues 
of which are complex. 

To verify this argument, we have computed the eigen- 
values of Pp^ for both horizontal and vertical flames. The 
results of the two cases are similar and thus we present only 
those of the vertical flame. We use the notation introduced 
by- Chong e/ al. 14 to classify the eigenvalues of P p i j . For the 
variable density case, the sum of the eigenvalues denoted by 
the invariant P is not necessarily zero, and thus the catego- 
ries of P>0 and P<0 exist. 

We find that, for the vertical flame, 85% of all mesh 
points belong to one of the eigencategories with all real ei- 
genvalues. Following Chong et al., 14 these categories with 
only real eigenvalues are labeled as 1 a (all negative eigen- 
values), lb (all positive eigenvalues), 6 a (two negative, one 
positive eigenvalue) and 6b (two positive and one negative 
eigenvalue). They contain, respectively, 6.6%, 7.2%, 35.72% 
and 34.9% of all the mesh points. The rest of the points 
(15%), are distributed among the categories 9a, 9b, 10a and 
10b which have one complex conjugate and one real eigen- 
value. (See Ref. 14 for details of the classification of eigen- 
values.) They contain respectively 4.9%, 2.8%, 3.1% and 
4.7% of all the points. 

Figure 12 shows the pdfs of the strain-enstrophy angle 
'P for each of these eigencategories. In this figure, the lines 
for categories with real eigenvalues which contain more 
points appear smoother. As was discussed earlier, in the F° 
band, the 'P values are close to 90°, and we see from Fig. 12 
that they belong to either la (i.e., ) or 6a (i.e., H ). 



FJG. 12. Strain-enstrophy angle Nk pdfs for different types of pressure 
Hessian eigenvalues. The labels of each group is taken from Chong et at. 
(1990). Top: Groups la, lb, 6a and 6b with only real eigenvalues. Bottom: 
Groups 9a, 9b, 10a and 10b with one conjugate pair complex eigenvalue. 
Vertical fiame. 


Figure 12 also shows that at the other extreme of the SP 
pdfs (i.e., 'l / ~0°) where vorticity dominates over strain, 
most points belong to the category lb (i.e., + + +). Between 
the two extremes of >P (i.e., ^P^O 0 and ^=90°), all the 
categories with complex eigenvalues are seen (Fig. 12) in 
addition to the group 6b which has all real eigenvalues. Re- 
call from Fig. 6(a) that mesh points with 'P~40° belong to 
the F + and F~ bands. 

We conclude that buoyancy produces complex P p i j ei- 
genvalues in the F + (and F~ for the vertical flame) band(s). 
On the other hand, the F° band (which has 'P es 90°) has 
mostly real P pi j eigenvalues indicating that buoyancy ef- 
fects are not considerable at the flame surface. 

V. CONCLUDING REMARKS 

In turbulent nonpremixed flames, buoyancy effects can 
be summarized as follows: Buoyancy generates strong vorti- , 
ces on both sides of the flame surface ( F sl ), and thus reduces 
the strain-dominance in the field. This reduction results in 
the vorticity alignment with the fi strain away from the flame 
surface. 

The strain-enstrophy angle 'P, can be used to determine 
the regions in which a or alignment trends are expected. 
The motivation for introducing 'P is based on the alignment 
equation [Eq. (1)] which we derived for the variable density 
case. This equation indicates that in regions where strain 
dominates over vorticity (or equivalently the baroclinic 
term), a alignment is expected. 

The alignment characteristics of the region near F„ (de- 
noted by F°) are not affected by buoyancy. This result is of 
importance to the scalar field. The peak scalar gradient (dis- 
sipation rate e F ) occurs in F° regardless of the presence of 
the buoyancy. 

Our results are consistent with the recent analysis of 
Gibbon and Heritage 12 and Galanti et al . 13 which suggested 
that the Navier-Stokes equations evolve toward an attracting 
solution. All the requirements for the solution of the Navier- 
Stokes equations to move to an attracting fixed point are met 
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in the flows wc considered in this paper. One drawback of 
the analysis of Refs. 12 and 13 is that it docs not distinguish 
between the a alignment and the /3 alignment. Also, the 
analysis 12,13 ignores the dynamics of regions with large 
strain/small vorticity. In fact, along a vorticity null line Ay 
and A P arc undefined. 
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